Current in narrow channels of anisotropic superconductors 
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We argue that in channels cut out of anisotropic single crystal superconductors and narrow on the 
scale of London penetration depth, the persistent current must cause the transverse phase difference 
provided the current does not point in any of the principal crystal directions. The difference is 
proportional to the current value and depends on the anisotropy parameter, on the current direction 
relative to the crystal, and on the transverse channel dimension. An experimental set up to measure 
the transverse phase is proposed. 
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In isotropic superconductors the supercurrent density 
is proportional to the gradient of the gauge invariant 
phase V<p = Vx + 27rA/0 o , where x 1S the phase of the 
order parameter ip = IV'le 1 *, O = nhc/\e\ is the flux 
quantum, and A is the vector potential: 
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M is the carrier mass. If the cross-section of a su- 
perconducting channel has small dimensions compared 
with the London penetration depth A, both j and Wip 
have the only nonzero components along the channel. 
For straight channels with the long dimension along x, 
dyip = d z ip = 0, i.e., the phase is constant in the trans- 
verse directions. 

This, however, is not the case for anisotropic supercon- 
ductors, where 
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Here, is the mass tensor, and summation is implied 
over repeated indices. It is convenient to normalize the 
masses: m ik = M ik / \M a M b M c ) 1 / 5 = M ik /M. Then the 
eigenvalues of rriik are related by m a mbm c = 1. In the 
uniaxial case which we consider for simplicity, m 2 a m c = 1; 
the masses then are expressed in terms of a single pa- 
rameter, the anisotropy ratio 7 2 = m c /m a : m a 
m c — y 4 / 3 . 

We invert Eq. (0) to obtain 
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which shows that the phase gradient Vy> and the current 
j are not parallel unless both of them point in a principal 
crystal direction. 

Consider a channel of a rectangular cross-section and 
denote by W and d its width in the y direction and the 
thickness in the z direction. To further simplify the prob- 
lem, we take z as one of principal crystal directions, say 



b. The axes c and a are then situated in the xy plane as 
shown in Fig. Denoting by the misalignment angle 
between the c axis and x, we readily obtain: 
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whereas m 



zz = m b = 7 
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FIG. 1: The current I in the superconducting channel of the 
width W (0 < y < W) and of the thickness d. The crystal axis 
b coincides with z. The cross denotes the Josephson junction. 

Let a small supercurrent / be fed into the channel along 
x. Then, Eq. (^|) yields the only non-zero transverse 
component of Vc/?: 



dyip 



M 



(5) 



2e?#| 2 yx Wd ■ 
Therefore, the side faces of the channel at y — 0,W 
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should possess the phase difference of 
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for any fixed x. 

In principle, this phase difference can be recorded by 
attaching a superconducting wire between the points 
{x, 0} and {x, W}. If the wire carries the Josephson junc- 
tion, the current in the wire is 
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where Iq is the maximum Josephson current. Thus, the 
current I w in the wire oscillates as a function of the driv- 
ing current / in the channel with the period 
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Near the critical temperature T c , the domain for which 
the above formulas are written, the equilibrium order pa- 
rameter %j) in clean materials is related to the carrier den- 
sity n e : \ip\ 2 ~ n e r, r = 1 — T/T c . This estimate holds 
provided / <C Id P where Id P « {ccf)Q / 1Q-k 2 X 2 ^)W d is the 
depairing current. 

To find out conditions under which the current oscilla- 
tions can be seen near T c , we require that both / and AI 
are small relative to Id p - Taking W ~ A(T) w Xo/y/r ~ 
Mc 2 / A-Ke 2 n e y/r and £ ~ £o/v^~j we obtain: 



m xy > 167t/k , 
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where k = Aq/£o- Since m 
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, c = 7 4 ' 3 , this inequality 
can be satisfied for strongly anisotropic materials with 
large Ginzburg-Landau parameter k. 

At low T's one can use the London expression for the 
current density instead of Eq. (0) : 
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where A = (A^Ae) 1 / 3 . We, therefore, can replace in the 
above formulas 2eh\ip\ 2 /M with c</>o/47r 2 A 2 . In particu- 
lar, we have 
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We then obtain a condition m xy 3> 8tt/k similar to (Q) 
in the whole temperature domain; if this condition is 
not satisfied, the window of driving currents / for the 
oscillations of I w to occur becomes narrow or disappears. 

Thus, the period A J is proportional to the channel 
thickness d and to t, the temperature distance from T c . 
Taking for an estimate realistic values of A = 2000 A, 
d = 1000 A, and 7 = 10, we obtain AI w 1 mA for low 
temperatures; the period shrinks on approaching T c . 



It is worth noting that for the effect to be observable, 
the wire contacts must be small relative to W . Otherwise 
the variation of the phase along the channel will destroy 
the quantum coherence. On the other hand, if the point 
contacts are fixed at different values of x, an additional 
phase difference 
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enters the argument of the sine in Eq. (Q); here Ax — 
\x2 and X\^2 are the contacts positions. This contri- 
bution changes the period AI of the Josephson current 
oscillations, but the very fact of periodic dependence of 
I w on / persists. 

If the driving current / changes with time, so does the 
phase difference Then, the junction in the wire is 

subject to a voltage 
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where lo is the driving frequency. For the parameters used 
above and u> ~ 10 9 Hz, the voltage may reach values of 
a few /iV; this estimate improves if thinner samples and 
higher temperatures are used: V oc 1 /rd. 

Writing the driving ac current as J = I a smujt, we 
readily obtain: 

Iw = lo siii(ip a sin ut) 

00 
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where J's are the Bessel functions and 
47r 2 A 2 m, 
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For small amplitudes I a of the driving current, one can 
keep only the term n = in the series (|l4|): I w = 
IofaSmuit, i.e., the junction response is linear. Since 
the voltage in this case is V — (4>of a /^c)oj cos tut, one 
can say that the junction is loaded with an inductance 
L = c4>q/2itIq. This simple interpretation does not hold 
if the condition (p a <C 1 is violated; the junction then 
should generate odd harmonics. 

We have considered a particular case of a misalign- 
ment between the supercurrent direction and the crystal. 
Of course, situations different from that of Fig. 1 can 
be readily treated. Our aim is to turn the community 
attention to the "transverse" phase difference which 
must accompany supercurrents as long as they are not 
parallel to principal crystal directions. We believe that 
effects related to the transverse phase in anisotropic 
materials and their complexity are not exhausted by 
simple examples we describe. In particular, we did not 
consider details of the electric field penetration into 
the channel (implicitly assuming the penetration length 



3 



> W). Given richness of time dependent Josephson 
phenomena, one may envisage a host of time dependent 
effects related to the transverse phase. It should also 
be stressed that phenomena we describe, happen on 
the scale of the London A, which might be large as 
compared to the coherence length £, the scale relevant 
to the well-studied phase slips. 
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